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ABSTRACT 

Criteria for gravitational collapse of expanding shells in rotating, shearing galaxy disks 
were determined using three-dimensional numerical simulations in the thin shell ap- 
proximation. The simulations were run over a grid of 7 independent variables, and the 
^\ ' resultant probabilities for triggering and unstable masses were determined as func- 

tions of 8 dimensionless parameters. When the ratio of the midplane gas density to 
the midplane total density is small, an expanding shell reaches the disk scale height 
' and vents to the halo before it collapses. When the Toomre instability parameter Q, or 

C*") . a similar shear parameter, Qa, are large, Coriolis forces and shear stall or reverse the 

collapse before the shell accumulates enough mass to be unstable. With large values 
of c b sh / (GL), for rms velocity dispersion c s h in the swept-up matter and shell-driving 
luminosity L, the pressure in the accumulated gas is too large to allow collapse during 
' the expansion time. Considering ~ 5000 models covering a wide range of parame- 

ter space, the common properties of shell collapse as a mechanism for triggered star 

/ 2\°- 5 

formation arc: (1) the time scale is ~ 4 I c s h/2iTGp [GL] ' ) for ambient midplane 

density p, (2) the total fragment mass is ~ 2 x 10 7 M Q , of which only a small fraction 
q ■ is likely to be molecular, (3) the triggering radius is ~ 2 times the scale height, and 

5-h ' the triggering probability is ~ 0.5 for large OB associations. Star formation triggered 

c/3 , by shell collapse should be most common in gas-rich galaxies, such as young galaxies 

C3 ■ or those with late Hubble types. 
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1 INTRODUCTION a sequence of primary and secondary star formation (see 

. . . . , review in Elmegreen 1998). 

Giant shells m nearby galaxies often contain young clusters „. r , . . . „ , . , , . . 

btar formation that is sequentially triggered like this 
on their periphery, suggesting that gravitational instabilities , , , , , , , . . „ . , . . 

should be able to continue lor an extended time as one gen- 
m swept-up gas lead to co- moving, self-gravitatmg clouds , , . . ,,,, ,, 

, . , _, . ' , r X, . eration leads to another m any remaining gas. Whether it 

that collapse into stars (McCray & Kalatos 1987; lenono- ,. . , „ ., , , nll , , , , , ,, 

J; , continues mdehmtely and falls a whole galaxy depends on the 

Tagle & Bodenheimcr 1988 , or that pre-existing clouds are , . . „ . „ „ , lri _„ „ . „ . 

° ' triggering efficiency (Mueller & Arnett 1976; Gerola & Sei- 



compressed to an unstable state as the shell passes, pro- 
ducing the same result (Woodward 1976; Klein et al. 1985; 
Dopita, Mathewson & Ford 1985; Boss 1995; Foster & Boss 
1996; Vanhala & Cameron 1998; Abraham, Balazs, & Kun 
2000; Yamaguchi 2001a). Sometimes old clusters are found 
inside the shells (e.g., Patel et al. 1998; Wilcots & Miller 



den 1978). If the mass of the second generation is larger than 
the mass of the first, then the activity should grow until it 
dominates all star formation. This may be the case in the so- 
lar neighborhood, where most star formation occurs in dense 
clusters (Lada, Strom, & Myers 1993; Carpenter 2000) that 

, , . , , , are located at the edges of high pressure regions, indicative 

1998; Steward et al. 2000; Stewart & Walter 2000; Yam- ,, . . , iV , t.-cj h ^ i t a \ e n ■ u 

, _._„,, s . , ' , , . of triggering (e.g., the Irmd Nebufa: Leffoch & Cermcnaro 

aguchi et al. 2001b), in which case the morphology suggests 20Q0 . Rogette: phelpg & ^ ^ ^.^ Reipurth; Rq _ 

driguez & Chini 1999; Ophiuchus: de Geus 1992; Sco Cen: 
Preibisch & Zinnecker 1999; Cepheus OB3 and Perseus OB2: 
* E-mail: bge@watson.ibm.com Sargent 1979; W3/4/5: Carpenter, Heyer,& Snell 2000; see 

t E-mail: palous@ig.cas.cz review in Elmegreen et al. 2000). Infrared sources around the 

+ E-mail: sona@ig.cas.cz supernova remnant G349.7+0.2 also look triggered (Reynoso 
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Figure 1. (bottom) The time at the beginning of the instability, 
t;, (open circles), and the time of fragmentation, tf (crosses), are 
shown versus the dynamical time in the midplane for a variety of 
cases with different parameters, (top) The time is multiplied by 
the dimcnsionlcss parameter (GL) ' 1 /c ; 5 for bubble luminosity 
L and rms velocity dispersion in the shell, c 3 h- Most collapse 
times follow the regular relations given by equations 6 and 7. 
Those with very high normalized tf had high external velocity 
dispersions and barely collapsed before they stalled. 



& Magnum 2000) , as do hot cores in the galactic center cloud 
Sgr B2 (Martin-Pintado et al. 1999). 

Here we investigate the conditions for shell trigger- 
ing. Numerical models of expanding, non-magnetic, self- 
gravitating shells are made, following the method of 
Ehlerova & Palous (1996). Other models are in Ehlerova 
et al. (1997), Efremov et al. (1999) and Ehlerova & Palous 
(2002). The models use the thin shell approximation pro- 
posed by Sedov (1959), in which a thin, 3-dimensional shell 
surrounding a hot medium is divided into a number of ele- 
ments and a system of equations for motion, mass and en- 
ergy is solved. The approximation is appropriate for a blast- 
wave propagating in the interstellar medium because the 
radius of the shell is much larger than its thickness during 
the majority of the evolution. The main advantage of this 
method over solutions of the complete hydrodynamical and 
Poisson equations in 3 dimensions is a much lower demand 
of CPU time. This gives us the possibility to explore an ex- 
tended grid of independent parameters. However, the adop- 
tion of this method restricts the spatial resolution, which is 
limited by the finite thickness of the compressed layer. 



Figure 2. The fraction of the models listed in Table 1 that frag- 
mented before they stalled is plotted versus the 8 dimcnsionlcss 
parameters. Curves with a strong monotonic variation indicate 
that the instability has a sensitive dependence on that parame- 
ter. The most important parameters are the Toomre parameter, 
Q, and a similar shear parameter, Qa, which determine the im- 
portance of rotation and shear, the density relative to the critical 
tidal density, p/ptidal, the normalized scale height, which is equiv- 
alent to the gas mass fraction in the disk, and the dimcnsionlcss 
rms speed in the shell, c s j,/ (GL) ' 2 . The figure in the top left in- 
dicates that once the time gets much larger than ~ 0.4 times the 
disk dynamical time, shell collapse becomes unlikely for typical 
parameters. 



2 MODELS 

The interstellar medium is modeled by a stratified, non- 
magnetic disk with a Gaussian density profile of dispersion 
H, midplane density p, velocity dispersion c ex t, local rota- 
tion curve V(r) oc R a , and local angular rate, fi. Three di- 
mensional expanding shells were driven into this medium by 
an energy injection rate L = E/T for energy E and time T, 
after which the shells coasted on residual internal pressure. 
Values of E = 10 54 ergs and T = 1.5 x 10 7 yr correspond to 
a medium-size OB association. 

The expected resistance to expansion by external mag- 
netic compressive forces can be modeled approximately by 
increasing c ex t to a value larger than the turbulent speed, 
perhaps by a factor of 2. The force from magnetic tension 
cannot be modeled this way; it tends to elongate a bubble 
or contain it within the disk (Tomisaka 1990, 1998; Ferrier 
et al. 1991). The implications of these and other magnetic 
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Figure 3. The number of models used to make Figure 2 is shown 
for each dimensionless parameter. 



Figure 4. The average total fragment mass at the time of signifi- 
cant collapse, tf , is shown versus the dimensionless parameters for 
each run that was unstable. The fragment mass does not depend 
much on any parameter if the shell goes unstable. 



forces on bubble instability are not known and will require 
more elaborate modeling. 

The effects of density and velocity irregularities in the 
swept-up gas and surrounding medium are not considered 
(see Silich et al. 1996). Gravitational instabilities develop 
faster in the denser parts of this medium, so the time scale 
for instabilities and triggering in a clumpy gas can be shorter 
than the time scale given here for a uniform gas. We ex- 
pect that the clump triggering time can still be estimated 
from our results using the derived density dependence for 
shell collapse at the higher density of the clump. The thin 
shell approximation breaks down for small dense clumps, 
however, because they erode by surface instabilities and 
get compressed in 3-dimensions when the blastwave passes 
(MacLow et al. 1994; Boss 1995; Xu & Stone 1995; Vanhala, 
& Cameron 1998; Fukuda & Hanawa 2000). Our model is 
not meant to address this situation. In a real shell, there 
may be localized triggering inside small clouds that are en- 
gulfed by the expansion (e.g., Yamaguchi et al. 2001a), in 
addition to large-scale triggering along the shell periphery 
where the swept-up gas collapses into new clouds. Here we 
consider the latter process. 

An expanding shell is approximated by an infinitesi- 
mally thin surface with a local expansion speed v, radius R, 
surface density E, and velocity dispersion, c sh . The time h 
at the beginning of the instability is given by (Elmegreen 
1994) 




I0 b 10 7 10" 

Total Unstable Fragment Mass 



Figure 5. A histogram of the total unstable fragment masses 
showing the tendency for all unstable shells to have a fragment 
mass of a few times 10 7 Mq, for a wide range of ambient densities, 
scale heights, and other parameters. 
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(1) 



For t > t b and uj(t) > 0, a fragmentation integral determines 
the time of significant collapse tf. 
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The Toomre Q parameter for epicyclic frequency k 
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Figure 6. The average ratio of the midplane shell radius to the 
scale height at the time of significant fragmentation versus the 
dimensionless parameters for each run that was unstable. This ra- 
dius is the distance of the most unstable part of the shell from the 
expansion center. The shell radius is always about 2 scale heights 
when stars begin to form. The dependence on c ex t / (GL) 02 comes 
from the theoretical expression for the expansion of a bubble. The 
dependence on Qa comes from the shear that stretches a bubble 
and increases R. 



If(tf) = / u{t')dt' 



(2) 



The simulation stops when v < c ex t everywhere in the 
galaxy symmetry plane. 

The galaxy and shell parameters were varied over the 
grid of values in Table 1 to search for trends in the begin- 
ning instability time, the collapse time, and the unstable 
mass and radius. The outcome of the expansion depends 
on more than one dimensionless variable because there are 
many competing processes that can prevent collapse. Galac- 
tic rotation introduces a Coriolis force that attempts to turn 
around the expansion of a bubble and limit its total radius. 
Galactic shear spreads out the shell, making an ellipse and 
concentrating the accumulation at the tips (Palous, Franco, 
& Tenorio-Tagle 1990). The finite scale height can lead to 
bubble blowout and depressurization. The external rms mo- 
tions lead to erosion when the bubble speed gets low. Ran- 
dom sonic and turbulent motions inside the swept-up matter 
stabilize it against collapse. 

To determine the effects of these processes, we consider 
various dimensionless parameters: 



SI (2 [1 + a]) ' and disk column density E = 2pH, 
Q tyGE 2irGpH 



(3) 



• An analogous shear parameter for Oort parameter A = 
-0.5RdQ/dR = 0.5D. (1 - a), 



»V3, 



tvGT, 



2-ivGpH ■ ^ 

The constant in this shear-related expression is designed to 
make Qa equal to Q for a flat rotation curve [a. = 0). 

• Another parameter involving gaseous self-gravity is the 
ratio of the midplane density to the density at the tidal 
limit. The tidal limiting density for a spherical cluster or 
cloud that co-rotates with the galaxy is pudai = 3AQ/ (nG). 

• The dimensionless scale height for total midplane den- 
sity ptotal, 

1/2 

(5) 



H (2nGp) 



1/2 



Cext 



Ptotal 



This expression uses H = c ext / (2nG ptotal) 1 ^ 2 , which is 
valid for an isothermal self-gravitating layer with p z (z) = 
ptotaisech 2 (z/H) . Our simulations use a Gaussian layer with 
p{z) oc exp (—z 2 /H 2 ), so the relation between E, c and H 
is not exactly the same. The difference between a Gaussian 
layer and a sech 2 layer is small for z < 2H, so we use the 
isothermal disk expressions in this paper. 

The dimensionless scale height is the square root of the 
ratio of the midplane gas density to the midplane total 
density and should be less than 1. In the simulations, H, 
c ext , and p are free parameters and some combinations give 
H (2-irGp) 1 ^ 2 I Cext larger than 1. 

• A dimensionless parameter that involves the internal 
rms speed and the shell luminosity is c s h/ (GL) 1 '' , and 
the same type of parameter with the external rms speed 
is Cext I (GL) 1 ^ 5 . The first of these luminosity parameters 
was studied by Whitworth & Francis (2002), who found 
c s h/ (GL) 1 ^ 3 < 0.13 as a condition for shell collapse. The 
second parameter was studied, in a slightly different form, 
by Palous & Ehlerova (2002), and Ehlerova & Palous (2001, 
2002). The ratio of the two rms speeds, c s h/c ex t, is also con- 
sidered here. 

The total number of parameter combinations from Ta- 
ble 1 is around 20,000, but most of these give unrealistic val- 
ues for the dimensionless parameters. We therefore limit the 
input parameter values to those which lead to 0.3 < Q < 15, 
H (2%Gp) 1/2 /cext < 3, and p/p tid ai < 1. 



3 RESULTS 

The times for the beginning of the instability, given by equa- 
tion jj], and the beginning of significant fragmentation, given 
by equation |2| are plotted in Figure 1 versus the midplane 
dynamical time, (2nGp)~ 1 ^ 2 . The range of parameters used 
for this plot differs from the range in Table 1 because here we 
include many different densities. Thus we consider Q — 0.5, 
1, 2, and 4 instead of the density explicitly, and we take 



Cext = 4 to 22 km s in steps of 2 km s 



and 



Csh 



0.5, 



1, and 3 km s . We also fix Q, = 25.9 km s kpc 
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Parameter Values 



Rotation curve slope, a -0.5, 0, 0.5, 0.8 

Angular rotation Rate, f2 5, 10, 25, 50, 100 km s kpc 

Bubble Energy, E 10 53 , 10 54 , 10 55 erg 

Injection Time, T 15 Myr 

Disk Scale Height, H 0.1, 0.2, 0.4, 0.8 kpc 

Ambient Density, p 0.1, 0.3, 1, 3, 10, 30, 100 H cm -3 

External rms velocity dispersion, c ex t 4, 8, 12, 20 km s — 1 

Internal shell velocity dispersion, c s ^ 0.3, 1, 3 km s — 1 



Table 1. Parameter Values 



in the Solar neighborhood. The parameter choices for H, a, 
and L are the same as in the table. 

Most of the scatter at the bottom of Figure 1 is from 
variations in shell rms speed, c s h, and bubble luminosity, L, 
which combine into the dimensionless parameter c\ h j (GL). 
The scatter is reduced at the top of the figure if we nor- 
malize time to (c a h/ [GL] ' 2 ) . The resulting normalized 
correlations are 

l0.2\ 0.5 



density is converted into a mass column density E = 2pH 
using the velocity dispersion to get the scale height, H — 

-1/2 j.. . a . 



c ext {2TtG P y L ", then tf ~ 4 (c sh /[GL]° 



^/(vrGE). 



t b ~ 1.5 



and 



ti 



shl [GL] 
2%Gp 



sh/ [GLf 
2-kGp 



(6) 



(7) 



The relative importance of the various competing pro- 
cesses is illustrated in Figure 2, which plots the fraction of 
the models in Table 1 that led to collapse versus each dimen- 
sionless parameter. The numbers of models used to deter- 
mine these fractions are shown in Figure 3 to illustrate the 
breadth of coverage of the parameter space and the likely 
uncertainty in the triggering fractions. The top left panel 
in Figure 2 has a lower triggering fraction than the other 
panels because the normalization for each value of the ab- 
scissa in the top left is the constant total number of runs, 
which is large. In the other panels, the normalization is the 
number of runs in each bin of the dimensionless parameter. 
We have to use the total in the top left because if there is 
no triggering, then tf and the abscissa are not defined. 

Figure 2 indicates that gravitational collapse in shear- 
ing, expanding shells is enhanced for large values of 
H (2-ivGp) 1 ^ 2 /cext, small values of Q and Qa, large values of 
p/pudai, and small values of c s h/ (GL) ' 2 . In the first case, 
the instability is unlikely for small relative gas densities be- 
cause the shell punches through the thin gas layer before it 
goes unstable. The instability is also unlikely for large Q and 
Qa because the corresponding large Coriolis force and shear 
slow the expansion, causing the shell to twist around and 
limit its accumulation of mass. Collapse is more likely for 
large p/pudai because the ambient medium is close to insta- 
bility anyway. Collapse is more likely for small c B h/ (GL) ' 2 
because then the shell is thin and dense for a given radius, 
so the pressure in the swept-up gas cannot easily resist its 
self-gravity. With c a h/ (GL) ' 2 < 0.02 for significant trigger- 
ing, the fragmentation time becomes tf < 0.5 (2nGp)~ 1 ^ 2 , 
which is less than half the dynamical time in the external 
medium. 

The density dependence in tf from equation (^) helps 
to explain the Q dependence in Figure 2. If the ambient 



This unstable time, multiplied by ft, is 4Q (c s h/[GL]°' 2 J ' . 

Thus large Q or large {c ah /[GL]°' 2 ^° ' 5 correspond to an un- 
stable time longer than an epicycle time, in which case the 
shell growth stops by Coriolis forces before it becomes unsta- 
ble. If c a h is small because of shell cooling, then the shell can 
collapse even if Q is large enough to make the disk otherwise 
stable. 

The total unstable mass in the shell at the fragmenta- 
tion time tf varies by only a small amount, a factor of ~ 10, 
for all of the unstable models. The average of this total mass 
scales mostly with the bubble luminosity, but does not cor- 
relate with any of the dimensionless parameters, as shown 
in Figure 4. Figure 5 plots a histogram of the number of 
shell models in regularly spaced logarithmic intervals of the 
total unstable mass. The shells need a sufficiently large mass 
to be unstable for typical expansion speeds. The total mass 
involved is ~ 2 x 10 7 Mq, although the molecular part that 
forms stars will be smaller. 

The shell radius at the fragmentation time, Rf, is shown 
relative to the scale height, H , in Figure 6. This radius is de- 
fined to be the distance from the expansion center, measured 
along the plane, of the most unstable part of the shell at the 
time of fragmentation. The ratio Rf/H depends mostly on 
the internal sound speed, aside from the usual dependence 
On Cext I (GL) ' 2 that comes into the theoretical expression 
for the bubble expansion (Weaver et al. 1977). There is a 
slight correlation between Rf/H and Qa, but not with Q, 
because shear stretches a bubble and gives it a larger average 
radius at the fragmentation time. The results suggest that 
the shell radius along the plane is ~ 2 to 3 scale heights when 
star formation begins. This distance corresponds to ~ 200 
to 300 pc for typical disks, in agreement with observations. 
The shell size should be ~ 2 times larger in dwarf galaxies 
because of their larger scale heights (e.g., Walter & Brinks 
1999). A shell is larger than this perpendicular to the plane 
because of its more rapid expansion into lower densities. 



4 CONCLUSIONS 

Sequential star formation triggered by the collapse of ex- 
panding shells requires essentially three conditions: that the 
midplane density of the gas be comparable to the stellar 
midplane density, that the Toomre parameter, Q, be rela- 
tively low, and that the velocity dispersion in the swept- 
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up gas be relatively low compared to (GL) ' 2 . The first of 
these conditions avoids blow out above and below the plane 
in the time it takes the shell to go unstable. The second 
avoids shell stalling by Coriolis forces, and the third makes 
the shell dense and strongly self-gravitating (Whitworth & 
Francis 2002). Other dimensionless parameters have a sim- 
ilar sensitivity: Qa and the inverse of p/pudal behave like 
Q. 

When and where are the dimensionless conditions for 
triggering fulfilled? In early type spiral galaxies where the 
relative gas density is low, Q, Qa and ptidai/p can be large 
(Caldwell et al. 1992) and H {2-nGp) 1 ^ 2 /c ext can be small. 
In this case, shell instabilities require very small values of 
c sh /(GL) - 2 ; i.e., either very large shell-driving luminosities 
or very low rms dispersions and temperatures in the swept- 
up gas. This is rather exceptional, although perhaps the 
triggering conditions are fulfilled in the vicinity of OB 78 
in M31. In late type spirals and gas rich dwarf galaxies, Q, 
Qa, and pudai/p are small and H (2nGp) 1/2 /c ext is large. 
Consequently, triggered star formation in expanding shells 
should be more frequent, as in the supergiant shells of the 
LMC (Yamaguchi et al. 2001a,b) and IC 2574 (Stewart & 
Walter 2000). 

The relative gas density is often high in starburst re- 
gions too, but there the absolute gas density is high as well, 
so the interstellar dynamical time can be very short. If it is 
shorter than the lifetime of an O-type star, which is typical 
for nuclear starburst regions, then bubble expansions will 
not be significantly aided by supernova explosions before 
the collapse time. This lowers L for a given OB association 
mass. Starburst regions can also have large c ex t, which limits 
shell growth and makes HII regions ineffective at moving gas 
around. As a result, spontaneous star formation processes, 
or those triggered by random supernovae, are more likely. 

Perhaps the best place for shell-triggered star formation 
was in the early Universe, when most galaxies were gas-rich 
and the dynamical time for an average galactocentric radius 
was much longer than an O-star lifetime. Galaxies at high 
redshift should contain many gaseous shells and triggering 
events, causing them to resemble sheared and massive ver- 
sions of today's irregular dwarf galaxies. 
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